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Abstract. A double Tavis-Cummings model (DTCM) is developed to simu- 
late the entanglement dynamics of realistic quantum information processing 
where two entangled atom-pairs AB and CD are distributed in such a way 
that atoms AC are embedded in a cavity a while BD are located in another 
remote cavity b. The evolutions of different types of initially shared entan- 
glement of atoms are studied under various initial states of cavity fields. The 
results obtained in the DTCM are compared with that obtained in the dou- 
ble Jaynes-Cummings model (DJCM) [J. Phys. B 40, S45 (2007)] and an 
interaction strength theory is proposed to explain the parameter domain in 
which the so-called entanglement sudden death occurs for both the DTCM 
and DJCM. 
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1 Introduction 



Entanglement is not only a key concept to distinguish between the quantum and the classical 
worlds, but has also been viewed as an indispensable resource to perform various intriguing 
global tasks in quantum computing and quantum information processing [1]. However, a 
notable characteristic of entanglement is its fragility in practical applications due to unavoid- 
able interaction with the environment. It is therefore of increasing importance to understand 
entanglement from its dynamical behaviors in realistic systems. As a rule for a global task, 
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entanglement should be shared between different remote parties who participate in the task. 
There are cases like teleportation [2], remote state preparation [3], etc., in which each particle 
of a multipartite entangled state is distributed to a separate location. There are also cases 
in which the entangled particles should be distributed so that each location contains several 
particles. For example, in quantum secret communication protocol between Alice and Bob 
[4], an ordered N Einstein-Podolky- Rosen (EPR) pairs are to be shared in such a way that 
Alice and Bob each holds one half of the pairs. That is, at Alice's location there are N 
particles which interact with one environment while the other partner-particles at Bob's 
location collectively interact with another environment. This scenario results in two inde- 
pendent local environments but each of them is common for one half of the A^ EPR pairs. A 
natural question arises as to how such kind of particle-environment interactions degrade the 
originally prepared global entanglement. This question is of fundamental interest because 
any quantum protocol depends essentially on the quality of the shared entanglement. As 
a first step to the problem, in this paper, we consider the case of A^ = 2 with two pairs of 
entangled two-level atoms AB and CD prepared in one of the two types of Bell-like states, 
namely, 

|^(0))/j = cos(a)|10)jj + sin(a)|01),j, (1) 

and 

\(p{0))u = cos(a)|ll),j + sin(«)|00),j, (2) 

where IJ G {AB,CD} and |0) (|1)) is the atomic ground (excited) state. 

For the simplest case of A^ = 1, i.e., either state ([T]) or state ([2]) is concerned for the 
initial state of a single atom-pair, the so-called double Jaynes-Cummings model (DJCM) 
[5-12] has been extensively adopted to study this problem because it yields exact analytical 
results. In the DJCM, each of two entangled atoms is embedded in an independent cavity 
and locally interacts with it. The results obtained within the DJCM for the initial empty 
cavities are that for any value of a state ([T]) loses its entanglement only at discrete time 
moments ti = {I + l/2)7[/g with / = 0,1,2,... and g the atom-cavity coupling constant, 
but for a certain domain of a state ([2]) may become separable at times smaller than ti 
and remains unentangled for some duration of time [6]. The latter phenomenon is referred 
to in the current literatures as entanglement sudden death (ESD) [13], which has been 
experimentally observed in [14,15]. An entangled state with ESD in evolution is less robust 



FIG. 1: Schematic representation of two entangled atom-pairs AB and CD which atoms A and 
C are located in cavity a but atoms B and D in another cavity b. 

than states without it, since ESD puts a limitation on the application time of entanglement. 
Therefore, studying ESD, especially conditions and parameter domains for its occurrence, 
is important from both theoretical and practical points of view. In Ref. [10] the DJCM 
is considered again and it is found that if the cavity fields are initially in Fock states with 
nonzero photon numbers then both atomic states |V^(0)) and \if{0)) would suffer from ESD 
for all values of a. The DJCM was also investigated from other perspectives and it was 
shown that the entanglement evolution of atoms is closely related to their energy variation 
[9] and there is a natural entanglement invariant demonstrating the entanglement transfer 
among all the system's degrees of freedom [7]. 

For the case of = 2 involving two pairs of entangled atoms, the situation would 
become more complex than that of = 1, because in each local environment there are two 
atoms simultaneously interacting with it. When there are many atoms interacting resonantly 
with a single-mode quantized radiation field of one and the same cavity the exact solution 
can be obtained by means of the so-called Tavis-Cummings model (TCM) [16]. Such a 
single TCM was used in Refs. [17] and [18] to study entanglement dynamics of two atoms 
that are initially prepared in a separable and entangled state, respectively. In this work 
we develop the so-called double Tavis-Cummings model (DTCM) including four two-level 
atoms A,B,C,D and two separate single-mode cavities a,b (see FIG. 1), which suffices for 
our purpose to study the entanglement dynamics for case of = 2. In the DTCM, atoms 
A (C) and B (D) are initially prepared either in state ([I]) or ([2]), but atoms A and C {B 
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and D) are located in cavity a (b) and interact with the cavity through the Tavis-Cummings 
Hamiltonian. We study the entanglement djTiamics of atom-pairs AB, CD, AC and BD by 
means of concurrence in dependence on the initial entanglement type of the atoms and on 
the initial state of cavity fields. We compare our results obtained in the DTCM with those 
obtained in the DJCM and present an interaction strength theory to explain the parameter 
domain in which the atom-pair exhibit ESD for both the DTCM and the DJCM. 

Our paper is organized as follows. In Sec. 2 we describe the DTCM and derive the exact 
analytical expression for the reduced density matrix of the atomic subsystem. Section 3 
presents detailed analysis of atomic entanglement dynamics when the initial atom-pairs are 
prepared either in state ([I]) or state (Ej) and the initial cavity fields are prepared either in 
the vacuum state, Fock state with a non-zero photon number or the thermal state. Finally, 
we conclude in Sec. 4. 

2 The double Tavis-Cummings model 

The total Hamiltonian of the system of four atoms A, B, C, D and two cavities a, b (see FIG. 
1) in the DTCM can be written as a sum of two isolated Tavis-Cummings Hamiltonians 

H = HACa + HsDh, (3) 

with 

HACa = ^(o-l + <^c) + ^(^'^0, + g {aaf + a+a,"), (4) 

and 

Hbd, = ^(^I. + ^L) + ^h-^h + g (bat + b^^-), (5) 

^ i=B,D 

where c^o {^^) is the frequency of the atom (cavity field mode), a (a~^) is the annihilation 
(creation) operator of the field in cavity a, b {b"^) is the annihilation (creation) operator of 
the field in cavity b, = (0| {cr^ = \0)-- (1|) is the rising (lowering) operator for the 
transition of atom i and g is the atom-cavity field coupling constant. Here, we are interested 
in the resonant case with ujq = u [16]. The initial cavity fields are assumed to be either in 
the vacuum state, the Fock state with a non-zero photon number or the thermal state. The 
general thermal field with its mean photon number n is a weighted mixture of Fock states 
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whose density operator pp can be represented as 

oo 

PF = ^Pn \n) {n\ , (6) 

n=0 

with \n) the Fock state of n photons and P„ is given by 

P = ^! (7) 

By virtue of the general thermal field defined above, through setting P„ = 5ni in Eq. ([2]), 
we can also study the vacuum state [l = 0) as well as any Fock states (Z > 0) of the fields. 
As for the initial states of atom-pairs AB and CD, we assume both of them to be either 
in state ([1]) or state ([2]). At t = the total state involving the four atoms and two cavities 
reads 

1 oo 
i,j,k,l=0 m,n=0 

where ao = sin a, ai = cos a and i'{j',k',l') = i{j,k,l) © 1 (with © an addition 
mod 2) for state ([1]), while i'{j',k',l') = i{j,k,l) for state ([2]). The evolution operator 
UACa(BDb){t) = exp{—iHACa{BDb)t) for the local interaction of atoms AC (BD) with cav- 
ity a (6) was derived exactly in Ref. [17]. At any time t > the state p(0) evolves into 
pit) = UACa(t)UBDbit)p{0)UXca{t)UBDbi^) which Can be represented as 

1 oo 

P(^) = E E aia.akaiP^P^ 

i,j,k,l=0 m,n=0 

UACait) \ik, (j7, m| U\caif) 

®UBDb{t) \i'k', n) BDbBDb 0''^'' '^l ^EDb- (9) 

Using the analytical expression of UACa{BDb)it) in [17] we have for UACa Vk,m) (similarly 
for UBDb \i'k\n)j^j^^) : 



UACa{t) \ik, m)^ca = E T)\i®p,k® q) h ~ ^-'^yP ~ ("^)'^). (10) 

p,q=0 



a 



where the functions Xik,pq{m, r) with t = gt are given in Appendix A for various possible 
i, k,p, q. These functions satisfy the normalization condition 

1 

Y |X,fc,p,(m,r)|2 = l (11) 

p,q=0 



6 



for any i, k, m and r. 

The reduced density matrix p^^'-'^{t) of the atomic subsystem can be obtained by tracing 
out p{t) over the cavity fields, i.e. 

^ABCD^^^ = Tr,,p(t) = ^ a,a,a,aiS^^c mACAc (j^l) ® 4d {f^) (12) 

where £^^y (K^)xyxy 0^1) ' with XFc = ACa or represents the map 

oo 

m,m'=0 
oo 1 

= H X! Pm^i-iyr— {-!)'' s,{-iyu-{-iyv 
m=0 r,s,u,v=0 

xXik,rs{m, r)X* _„^(m, r) |i © r, /c © s)^^^^ (j © m, / © t;| . (13) 

The exphcit expressions of Sxy i\'^k)xYXY 0^1) given in Appendix B for various possible 
i, k, j, I. 

3 Atomic entanglement dynamics 

With the formulae derived in the previous section we are now in the position to analyze 
the entanglement dynamics of any atom-pair. By using Eq. (|T2l) we can readily get the 
reduced density matrix of any pair of atoms by tracing out p^^'^^(t) over the degrees of 
freedom of the remaining atoms. In two-qubit domains, there exist a number of good 
measures of entanglement such as concurrence [19] and negativity [20]. Although the various 
entanglement measures may be somewhat different quantitatively [6], they are qualitatively 
equivalent to each other in the sense that all of them are equal to zero for unentangled 
states. Here we adopt Wootters' concurrence [19] because of its convenience in definition, 
normalization and calculation. The concurrence C for any (reduced) density matrix p of 
two qubits is defined as 

C(p) =max{0,yA^- y^- ^aI}, (14) 

where Aj (Ai > A2 > A3 > A4) are the eigenvalues of the matrix ( = p{ay © ay)p*{ay © a^), 
with ffy a Pauli matrix and p* the complex conjugation of p in the standard basis. For 
separate states C{p) = 0, whereas for maximally entangled states C{p) = 1. In particular. 
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if p is of the X-form [21], 



P 



IJ 



(15) 



qH qH 
qH qH 
\Qii Q qH I 

where g^j. are real positive and q{{ = (ftV) generally complex, then the concurrence 
( IT^ simplifies to 



C^-^ = 2max{0, \gii\ - Jgi(gli, \gii\ - x/gMgii}. 



JJ\ 



IJ jj \JJ\ 



(16) 



Since both states ([T]) and ([2]) of the atoms take on and preserve the X-form in their evolution, 
Eq. f|T6l) is very useful throughout this work. 



3.1 1^(0)) type initial state for atom-pairs AB and CD 

We first consider the case when both the atom-pairs AB and CD are initially prepared in 
state ([1]). In accordance with Eq. (fT2l) the reduced density matrix of the atomic subsystem 
at any time t is 

1 

pjBCDf^^^ = ^ a^ajakaiElc (K, k)ACAc 0' ^1) ® 4d (K © 1, A; © ^)bdbd (j © 1, / © 1|) , 

(17) 

which can be evaluated straightforwardly via the map ( fT3l) . Then the reduced density matri- 
ces of interest are p^^(t) =TrcDP/ ""^^(t), pF''(t) =TrABP^^^^^(t), pf{t) =TTBDpf''^''{t) 
and pf^(t) =TrAcP/^^^(t). All of pf^{t), pf'^{t), pf^{t) and pf^(t) have the X-form so 
the corresponding concurrences are determined by Eq. (fT6|) . In the following we study the 
time dependence of these concurrences for the fields in cavities a and h being initially in the 
vacuum state, the Fock state with a non-zero photon number or the general thermal state, 
respectively. 

In FIG. 2 we plot Cf^ (the same for Cf^ due to symmetry) as functions of rescaled 
time gt and a for the initial cavity fields being in the vacuum state. From FIG. 2 it is 
transparent that Cf^ vanishes after a finite time of evolution and remains zero for some 
period of time before increasing again. This dynamics holds in the whole range of a. A 
comparison between the DTCM and the DJCM [6] for the same initial preparation of the 
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FIG. 2: The concurrence Cab = Cf (t) as functions of rescaled time gt and a for initially both 
cavity fields are in the vacuum state and both atom-pairs AB and CD are in the |^(0))([I]) type 
state in the DTCM. 
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FIG. 3: The concurrence Cab = Cj{t) as a function of gt for various values of a for the same 
initial preparation of cavities and atoms as in Fig. 2 in (a) the DTCM and (b) the DJCM. 

cavities and atoms is shown in FIG. 3. Within the first cycle of evolution, in the DJCM 
(see FIG. 3b) Cf^ vanishes at the moment to = '^/{'^q) and grows up again right after to? 
while in the DTCM (see FIG. 3a) Cf^ = at a time shorter than and remains so for 
some time before reviving. This indicates that for one and the same empty cavity fields, 
|'0(O)) type initial state of atoms does not undergo ESD in the DJCM but it does in the 
DTCM. Therefore, the atomic entanglement dynamics is model-dependent apart from the 
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entanglement type itself. The physical interpretation behind such a clear distinction in the 
dynamical behaviors between the two models can be thought of as follows. If the cavities 
are empty, atoms in the ground state |0) remain unchanged and only atoms in the excited 
state |1) can interact with the cavity fields. Denoting by N\i) the number of atoms that 
may be populated in state |1) , the system-environment interaction can be classified into 
two regimes, "strong" and "weak" interaction regimes, depending on relative magnitudes of 
P> and P<, where P> (P<) is the probability that N\i) > {N\i) < N^) with A^^^ the number 
of cavities. In the DTCM considered here and the DJCM considered in [6,7] it is clear that 
Nc = 2. We define the following convention: the strong interaction regime corresponds to 
P> > P<, while P> < P< implies the weak interaction regime. In the DJCM the total 
system state of two atoms A, B and two cavities a, 6 at t = reads 

I^(0))ab I00)a6 = cosa|10)A,|00)Bfe + sin a|00)Aa| 10)66, (18) 

whereas in the DTCM the total system state of four atoms A, B, C, D and two cavities a, h 
at it = reads 

\m)AB \m)cD |00),, = cos^ a \mj,ca ^301. + COS « siu « 1 100)^^, 1010)^^, 

+ sinacosa |010)^^, llOO)^^, + sin^ a |000)^^, |110)W-19) 

From Eq. (fTS!) it follows that there is always only one atom (namely, either atom A in the 
first term or atom B in the second term) being in state |1) regardless of the value of a. That 
is, P< = 1 < P> = 0, resulting in the weak interaction regime in the DJCM for the whole 
range of a. However, what is followed from Eq. ( fT9|) is that for any value of a there are 
always two atoms (namely, either atoms A and C in the first term or atoms A and D in the 
second term or atoms C and B in the third term or atoms B and D in the fourth term) 
being in state |1) . That is, P> = 1 > P< = 0, resulting in the strong interaction regime in 
the DTCM regardless of the value of a. Therefore, it can be said that, when the cavities 
are initially prepared in the vacuum state, |V'(0)) type initial state of atoms exhibits ESD 
in the strong interaction regime (i.e., in the DTCM) but it does not in the weak interaction 
regime (i.e., in the DJCM), independent of the parameter a. 

The case when the initial cavity fields are in a Fock state with a certain nonzero photon 
number is illustrated in FIG. 4. In this case not only atoms in state |1) but also atoms in 
state |0) , i.e., all the present atoms, can interact with the cavity fields so that the interaction 
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FIG. 4: The concurrence Cab = Cf (t) as a function of gt for q = 7r/4 for initially the cavity 
fields are in different Fock states \mn)^jj and atom-pairs AB and CD are in the |V'(0)) ^ type 
state in the DTCM. 
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FIG. 5: The concurrence Cab = Cj{t) as a function of gt for a = it/A for initially cavity fields 
are in the thermal state with different mean photon numbers m, n and atom-pairs AB and CD 
are in the |V'(0)) ([I]) type state in the DTCM. 

regime is always strong resulting in ESD for whatever values of a. A remarkable feature is 
that Cf^ decays quicker and reaches zero in a shorter time for a larger initial number of 
photons in the cavities. The underlying physics for that feature is the intensification of the 
system-environment effective interaction with the increase of photon number contained in 
the cavities. 
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FIG. 6: The concurrence Cbd = Cf (t) as functions of gt and a for initially the cavity fields are 
in the Fock state 111)^6 and both atom-pairs AB and CD are in the state ([I]) in the DTCM. 




gt 

FIG. 7: The concurrence Cbd = Cf^{t) as a function of gt for various values of a with the same 
initial preparation of cavity fields and atom-pairs as in Fig. 6 in the DTCM. 

Figure 5 plots the evolution of Cf^ for the cavity fields being initially in the thermal 
state. The entanglement dynamics looks chaotic due to the nature of the thermal fields. As 
can be seen from FIG. 5, the larger the mean photon number (corresponding to the higher 
temperature) the shorter the death time of Cf^ and the longer its revival time. 

At this point let us study the dynamics of the two atoms that are located in one and the 
same cavity. These are atoms A and C in cavity a and atoms B and D in cavity b. Such 
atoms in the same cavity are absolutely uncorrelated at the beginning and also there are no 
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FIG. 8: The concurrence Cbd = Cf (t) as functions of gt and a for initially the cavity fields are 
in the thermal state with the mean photon numbers m = n = 1 and both atom-pairs AB and CD 
are in the state ^ in the DTCM. 

direct interactions between them during the entire course of evolution, in accordance with the 
problem Hamiltonians (jlj) and ([5]). However, an effective (indirect) atom-atom interaction 
is induced for t > thanks to the coupling of both atoms with a common environment. 
Such an effective atom-atom interaction could nontrivially affect their global behaviors. In 
fact, as investigated in Ref. [17], if the initial atoms are prepared either in state |01) or 
1 10) (1 11)), then they always get entangled with each other (remain unentangled) regardless 
of the nature of the cavity fields. But, if the atomic initial state is |00), then the field in 
the vacuum state leaves the atoms unentangled and the field in a Fock state with a non- 
zero photon number or thermal state can entangle them. Here, in the DTCM, at variance 
with the situation considered in Ref. [17], at t = the atoms in a cavity, though being 
independent of each other, are entangled with other atoms in another cavity. That is, we have 
at if: = in cavity a (b) a mixed state pf^^O) =TTBDpf^'~^^{0) = Z]ij=o'^f'^| K'i)yiCyic (^'^1 
(pf^(O) =TTAcpf''^''{0) = E-,=oafa| K © 1,J © ^)bdbd © © instead of a pure 
state as in Ref. [17]. Figure 6 plots the concurrence Cf^ as functions of gt and a with 
the initial fields in both cavities containing just one photon. This figure shows that the 
entanglement dynamics of the atoms is sensitive to a, as it should be. For example, in 
the region of a G [0,0.297r] atoms B and D can get entangled, but for a around 7r/2 no 
entanglement is generated through the whole evolution. These results are in full agreement 
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with those reported in Ref. [17] where a = (i.e., pf'^(O) = [00)^^^^ (00|) and a = n/2 
(i.e., pf^(O) = \^^)bdbd concerned. To get more insight into the effect of a on 

atomic entanglement generation we show in FIG. 7 a 2D plot of Cf^ as a function of gt with 
the initial cavity fields in the Fock states |1, l)^j, for various values of a. When a = (i.e., 
(0) = |00)^^^^ (00|), the entanglement of B and D emerges immediately from t = 0. 
Nevertheless, when a > the atoms remain unentangled for some initial period of time and 
suddenly become entangled at some later time. The larger the value of a the longer the delay 
time of entanglement generation. Such phenomena of delayed entanglement during the time 
evolution can be called "entanglement sudden birth" (ESB) [22] . The effect of thermal fields 
on inducing entanglement between atoms B and D is drawn in FIG. 8 with the cavity mean 
photon numbers m = n = 1, which agrees well with the result in Ref. [17] for a = 0. Since 
the thermal state is a weighted mixture of Fock states (see Eq. ([6])), it is a chaotic state with 
minimum information and so its effect is generally irregular. In comparison with the case of 
"corresponding" Fock states |1, 1)^^, one sees that the region of a allowing entanglement of 
atoms is much shrunk and the amount of generated entanglement is very small. The plots 
of C/^ can be obtained from those of Cf^ by making a change a — > a + 7r/2. 

3.2 1(^(0)) type initial state for atom-pairs AB and CD 

We next consider the case when both atom-pairs AB and CD are initially prepared in 
state ([2]). In accordance with Eq. (IT^ the reduced density matrix of the atomic subsystem 
at any time t is 

i,j,k,l=0 

In FIG. 9 we plot Cfj^ (the same for Cf/^ due to symmetry) versus gt and a for the initial 
empty cavity fields. It is visual from this figure that ESD occurs but not in the whole range 
of a, in clear contrast with the case shown in FIG. 2 when both the atom-pairs AB and CD 
are initially prepared in state ([T]). To derive the constraint on a that triggers ESD let us 
look at the total system state at t = : 

\m)AB b(0))cD |00)., = cos^ a im^ca \m BDb + COS siu 1 100)^^, |100)^^, 

+ sin a cos a |010)^^, |010)^^, + sin^ a |000)^^„ |000)^^,(21) 
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Obviously, the probability that all the four atoms are in state |1) is cos^a, the probability 
that only two atoms (namely, either atoms A and B or atoms C and D) are in state |1) 
is 2cos^asin^a and the probability that none of the atoms are in state |1) (i.e., all the 
atoms are in state |0)) is sin^a. That is, P> = cos^a + 2 cos^ a sin^ a and P< = sin^a. As 
mentioned in the previous subsection, the condition for the occurrence of ESD is that the 
interaction regime is strong, i.e., P> > P<- So, the values of a for which ESD occurs should 
satisfy the constraint 

sin^cK^. (22) 
v2 

Noticeably, this constraint is not coincident with that one in the DJCM for which the initial 
total system state reads 

\^{^))aB I00)afe = COS«|10)Aa|10)m + siu «|00)Aa|00)Bfe. (23) 

As followed from Eq. ( !23l) . the probability that the two atoms are in state |1) is cos^ a 
and the probability that none of the atoms are in state |1) is sin^ a. That is, P> = cos^ a, 
P< = sin^ a and thus the values of a, for which the system-environment interaction regime 
is strong (i.e., ESD occurs) in the DJCM, satisfy the constraint 

sin^cK^. (24) 

The constraints (122|) and fl2^ imply that the a-parameter domain in which the atoms suffer 
from ESD is wider in the DTCM than in the DJCM. 

The case for the initial cavity fields being in a Fock state |ll)Qf, is plotted in FIG. 10. 
A remarkable feature as compared with the vacuum fields case in FIG. 9 is that here ESD 
occurs in the whole range of a. Again, the physical reason for this is that in the presence of 
initial photons all the atoms are in interaction with the cavity fields (i.e., not only atoms in 
state |1) but also those in state |0) interact with the cavity fields). 

In FIG. 11 we plot CfP as a function of gt for the initial fields in a thermal state with 
different mean photon numbers for a given value of a. Comparing FIG. 11 with FIG. 5 
signals that with relatively small mean photon numbers (e.g., m = n = 0.1) the signature 
of ESD is less pronounced for the case when the initial atoms are prepared in state ([2]) than 
in state ([H). 

The entanglement generation dynamics of the atomic pairs AC and BD is similar to the 
case considered in the preceding subsection and thus will not be iterated here. 
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FIG. 9: The concurrence Cab = Cfj {t) as functions of gt and a for initially both cavity fields 
are in the vacuum state and both atom-pairs AB and CD are in the 1^(0)) ([2]) type state in the 
DTCM. 




FIG. 10: The concurrence Cab = Cfj (t) as functions of gt and a for initially the cavity fields are 
in the Fock state lll)^^ and both atom-pairs AB and CD are in the |v(0)) ^ type state in the 
DTCM. 

4 Conclusion 

In conclusion, we have, by means of concurrence, studied the entanglement dynamics of 
the DTCM motivated by certain realistic quantum information processing. The system is 
composed of four two-level atoms A, B, C, D and two spatially separated single-mode cavities 
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FIG. 11: The concurrence Cab = Cfj {t) as a function of gt for a = tt/A for initially the cavity 
fields are in the thermal state with different mean photon numbers m, n and both atom-pairs AB 
and CD are in the \ip{0)) ^ type state in the DTCM. 

a,b. Initially, atom-pairs AB and CD are prepared either in Bell-like state \ip{0)) ([1]) or 
Ifi^)) 0, while both cavities are prepared either in the vacuum state, the Fock state with 
non-zero photon numbers or the thermal sate. Independent atoms A, C {B, D) that belong 
to different entangled atom-pairs are embedded in one and the same cavity a (b) and interact 
with it through the Tavis-Cummings Hamiltonian. 

For the vacuum fields the \ip{0)) type initial state of atom-pairs AB and CD displays 
ESD for the whole value range of the parameter a which represents the initial entanglement 
degree of AB and CD. This result is in sharp contrast with the DJCM for which ESD 
does not occur at all for whatever values of a [6,7]. As for the \<f{0)) type initial state of 
atom-pairs AB and CD, ESD only occur for the value of a such that sin^ a < l/\/2, which 
is wider than that in the DJCM where ESD occurs just for a such that sin^ a < 1/2 [6,7]. 
Physically, these results (i.e., the domain of a for which ESD occurs) in both the DTCM 
and DJCM can be explained via the interaction strength theory according to which ESD 
occurs (does not occur) in the strong (weak) system-environment interaction regime. The 
interaction regime is identified by the number of atoms that can have interaction with the 
cavities, which is determined by the relative magnitudes of P> and P< defined in subsection 
3.1. Remarkably, the interaction strength theory turns out to apply also for the so-called 
triple Jaynes-Cummings model [23] for GHZ-like atomic states as well as for the case of 
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multiple dissipative environments with multiqubit GHZ-like atomic states [24,25]. 

We have shown that the non-vacuum environments of cavities have great effects on the 
appearance of ESD for atoms. That is, when the cavity fields are initially in the Fock 
state with a non-zero photon number or the general thermal state, ESD always happens for 
atom-pairs AB and CD regardless of the entanglement type they are prepared. Moreover, 
the more photon number in the Fock state or the greater the mean photon number in the 
thermal state the quicker the entanglement decay rate, i.e., the sooner the time of ESD 
occurrence. In terms of the interaction strength theory, these properties are explained by 
the physical fact that in the presence of nonzero (mean) photon number the interaction 
regime is always strong because all the atoms (i.e., not only those in the excited state as in 
the case of empty cavities) can interact with the fields. Thus, the actual system-environment 
interaction strength is now identified by the number of excitation which in these cases is 
proportional to the total number of both atoms and photons. 

We have also studied creation of entanglement between initially uncorrelated atoms A 
and C in cavity a {B and D in cavity b). Compared to the case of a = considered in 
Ref. [17] here we showed that for a 7^ there appears the so-called entanglement sudden 
birth, i.e., the formation of atomic entanglement does not take place at once as the system 
evolves but emerges suddenly at some delayed time, which is dependent on the value of a. 
The DTCM presented in this work could be extended to the general multiple case where two 
groups of multipartite entangled atoms are distributed in such a way that every two atoms 
from different group are located in the same environment. In this way, we can study not only 
the pairwise entanglement of atoms between any two nodes (cavities or local environments) 
via concurrence but also the entanglement of any atomic bipartition by means of negativity. 
These studies can reveal the degraded properties of various multipartite entangled state and 
thus be useful for the large-scale quantum information processing. 

Z.X.M. and Y.J.X. are supported by National Natural Science Foundation of China under 
Grant No. 10774088. N.B.A. acknowledges support from a NAFOSTED project No. NCCB- 
2009 and from the KIAS Scholar program. 



18 

APPENDIX A: THE EXPLICIT EXPRESSIONS OF Xifcpg(m,r) 

The functions Xjfc pq(m, r) appearing in Eq. (fTOj) for all possible i, k,p, q are given by 

TTl -\- 1 / 

Xn,ooKr) = ^— -^[cos(v/2(2m + 3)r) - 1] + 1, (Al) 



m + 1 



A:ii,io(m,r) = Xii,oi (m,r) = -i^ sin(^2(2m + 3)r), (A2) 



m + l)(m + 2) 



""^ = 2m + 3 ^[cos(V2(2^ + 3)r) - 1], (A3) 



m 



Xoi,w{m,T) =Xio,oi(m,r) = -z^^^^— ^ sin(^2(2m + l)r), (A4) 

Xoi^oo{m,T) = Xio,oo("^,r) = i[cos(y^2(2m + l)r) + 1], (A5) 

^01,11 (m, r) = Xio,ii(m, r) = ^[cos(^2(2m + l)r) - 1], (A6) 

/ 'jYi -\- 1 / 

A:oi,oi("^, 1") = A:io,io("^, r) = -i\ — — sin( J2(2m + l)r), (A7) 

y 2{2m +1) * 



^00,11 ("^,r) 



\ mim — 1) / 

^ [cos(^2(2m-l)r) - 1], (A8) 



Aroo,oi("^,r) = Xoo,io("^,r) = -i J^^^— — sin(y^2(2m - l)r) (A9) 



and 



^oo,oo(rn, r) = _ _-[cos(Y'2(2m - l)r) - 1] + 1. (AlO) 



APPENDIX B: THE EXPLICIT EXPRESSIONS OF S^^y {\ik)^yxy (j7|) 

The expressions of the map £xy (K^)xyxy 0^1) ' with XYc = ACa or BDb, appearing 
in Eq. (fT3!) for all possible i, k,j, I are given by 



S'xvimxYXYm) = J2 P^[\Xoo,nim,T)\'\ll)xYXY{n\ 

m=0 

+ |Xoo,ioKr)niO)xyxy(10| 

+Xoo,io("^, ^)^oo,oi ("^^ ^) 1 10)xyxy (01 1 
+Xoo,oi('Ti, r)Xoo io(m, r)|01)xyxy(10| 
+ |Xoo,oi(^,r)n01) XYXY (Oil 

+ |Xoo,oo(m,r)noO)xyxy(00|l , (Bl) 



(oiir 

oo 

XYXY ((01| + (10| 

m=0 

+Xoi,u{^, T)X*o,oo{m, T)|10)xyxy(00| 
+Xoi,oo(m, T)X*o_oo(m, T)|01)xyxy (00|] , 

^xy(|10)xyxy(00|) = (|00)xyxy(10|)* 

= E [^io,oi(m,T)X*o,io(m,T)|ll)xyxy(10| 

m=0 

+Xio,oi (m, r)X*o 01 (m, r) 1 1 1 )xyxy (01 1 
+^10,00 (m, r)X*o oo("^, r)|10)xyxy (00| 
+Xio,ii(m,r)X*o,oo(^,T)|01)xyxy(00|" 

£^^(|ll)xyxy(00|) = (|00)xyxy(ll|)* 

oo 

xyxy{00\, 

m=0 

CXD 

^xy(|0lWy(01|) = E^m[l^oi,ioKT)nil)xyxy(ll| 

m=0 

+Xoi,ii (m, r)X*i^oo(m, r) 1 10)xyxy (01 1 
+ |Xoi,ii(m,T)niO)xyxy(00| 
+Xoi^oo{m, t)X*i ^i(m, r) |01)xyxy (10| 
+ |Xoi,oo(m,T)noi) xyxy (01 1 
+ |Xoi,oi(m,r)n00)xyxy(00| 

^XY 

(|10)xyxy(01|) = £:^^ (|01)xyxy(10|)* 

oo 

= ^P4[|Xio,oiKr)nil)xyxy(ll| 

m=0 

+Xio,oo(m, T)X*i,n("^, r) |10)xyxy (10| 
+|Xio,ooKT)niO)xyxy(01| 
+ |Xio,ii(m,r)noi)xyxy(10| 
+Xio,ii(m, T)X*i_oo(m, T)|01)xyxy(01| 
+|Xio,io(m,T)nOO)xyxy(00|" 
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(|11> XYXY 



(Oil) = E^^y{\oi)xYXY{n\r 

oo 

= E [Xn,ooKr)X*i^,i(m,r)|ll)xyxy(10| 

m=0 

+Xn,oo(m, t)X*i no(m, T)|ll)xyxy (01| 



+Xn,oi {m, t)X*i^oi (m, r) 1 10)xyxy (00| 

+Xn,io(m, r)X*,_o,(m, r)|01)xyxy (00| 



(B7) 



^XY {\10)XYXY 



(10|) = $:P^[|Xio,oi(m,T)nil)xyxy(ll| 

m=0 

+ |Xio,ooKr)niO)xyxy(10| 
+Xio,oo(m, t)X*o,ii(^, T)|10)xyxy(01| 
+Xio,ii(m, t)Xio oo(m, r) |01)xyxy (10| 
+ |Xio,ii(m,r)n01) xyxy (01 1 + 
|^io,ioKT)nOO)xyxy(00| 



(B8) 



^XY (111) XYXY 



and 



^xy (1 11) xyxy 



(10|) = £:i^(|10)xyxy(ll|)* 

oo 

= ^ [Xn,oo(m,T)X*o,ooKT)|ll)xyxy(10| 

m=0 

+Xii,oo(m, t)X*o 11 (m, r) | ll)xyxy (01 1 

+Xn,oi (m, r)X*o_io (m, r) 1 10)xyxy (00| 

+ ^ii,io(m, r)X*o io(m, r) |01)xyxy (OoH (B9) 



(111) = ^P^[|Xn,ooKT)nil)™y(ll| 

m=0 

+ 1^11,01 (m,T)|2 (|10) + 101))^^^^ ((10| + (Oil) 
+|^ii,ii(m,T)nOO)xyxy(00|l . (BIO) 
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